Evidence of columnar order in the fully frustrated transverse field Ising model on the square lattice 
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Using extensive classical and quantum Monte Carlo simulations, we investigate the ground-state phase dia- 
gram of the fully frustrated transverse field Ising model on the square lattice. We show that pure columnar order 
develops in the low-field phase above a surprisingly large length scale, below which an effective U(l) symmetry 
is present. The same conclusion applies to the Quantum Dimer Model with purely kinetic energy, to which the 
model reduces in the zero-field limit, as well as to the stacked classical version of the model. By contrast, the 
2D classical version of the model is shown to develop plaquette order. Semiclassical arguments show that the 
transition from plaquette to columnar order is a consequence of quantum fluctuations. 



The Ising model in a transverse field, introduced by De 
Gennes in the context of ferroelectric systems [1] and defined 
by the Hamiltonian 

{'7> j 

where cr" are Pauli matrices, F the transverse field and 7^ 
Ising coupling constants, has imposed itself as one of the most 
useful minimal models to study quantum fluctuations [2]. In 
particular, the simple ferromagnetic version of the model un- 
dergoes, when the transverse field is increased, a phase tran- 
sition to a paramagnetic phase that embodies several of the 
basic aspects of quantum phase transitions [3]. Frustrated ver- 
sions of the model, i.e. versions in which it is impossible to 
minimize the energy for all bonds simultaneously, have been 
studied over the years [4], with in general very rich phase di- 
agrams [2]. 

Recently, special attention has been payed to the subclass 
of Fully Frustrated Transverse Field Ising Models (FFTFIM) 
[5, 6]. In these models, the couplings are limited to nearest 
neighbors and have equal magnitudes, but their signs are such 
that it is impossible to satisfy simultaneously all bonds around 
any elementary plaquette, resulting in a very strong ground 
state degeneracy. The AF Ising model on the triangular lattice 
is a typical example. On the square lattice, to which this paper 
is devoted, this can be achieved by choosing antiferromag- 
netic (positive) couplings (7,y - J) on every second vertical 
line, and ferromagnetic (negative) couplings (7,^ - - J) for all 
other bonds, resulting into an Ising model often referred to as 
Villain's odd model [7]. 

The interest in this family of models has raised because 
of their connection to the Quantum Dimer Model (QDM) of 
Rokhsar and Kivelson [8]. First introduced in the context of 
the Resonating Valence Bond (RVB) theory of high T^. cuprate 
superconductors, the QDM has become one of the paradig- 
matic model in the field of quantum spin liquids after Moess- 
ner and Sondhi have shown that, on the triangular lattice, it has 
an RVB phase [9]. On the square lattice, the QDM is defined 



by 

^^QDM = Xi["KI-><n|+h.c.) + 

v(|n><n:| + |n><ii|)], (2) 

(in terms of hard core dimer objects II) where t and V are 
the amplitudes of kinetic and potential terms, and the sum runs 
over all elementary square plaque ttes. In the limit V/J — » 0, 
and up to a trivial degeneracy associated to a global spin flip, 
the FFTFIM maps onto the purely kinetic QDM with V/t - 
[5,6]. 

Over the last few years, the properties of the QDM on sev- 
eral lattices have been investigated in great details, and in most 
cases a consensus has been reached regarding its phase dia- 
gram as a function of V/t [10], with one noticeable excep- 
tion: on the square lattice, while it is clear that the QDM 
has a columnar phase for V — » -oo and a staggered phase 
for V/t > 1, the situation between these two limits is still 
controversial. It was suggested on the basis of exact diag- 
onalizations (ED) that the system has columnar order up to 
V/t = 1 [11], but subsequent investigations have challenged 
this picture. ED on larger systems have been interpreted as 
evidence of plaquette order in the range (V/Oe < V/t < 1 
with (y/r)e - -0.2 [12]. The same conclusion has been 
reached with Quantum Monte Carlo (QMC), with a critical 
ratio (V/t)c - 0.6 significantly closer to 1 though [13]. Fi- 
nally, a new scenario has been put forward on the basis of ED 
and Green's function QMC according to which a mixed phase 
is present between V/t ^ and V/t ^ 0.8 - 1, in which all 
symmetries broken in the columnar and plaquette phases are 
simultaneously broken [14]. This scenario has been supported 
by QMC on a related hard-core boson model [15]. 

A major problem in resolving the issue is that there is no 
simple limit where the plaquette phase is known to exist for 
sure in the QDM on the square lattice. In this Letter, we in- 
vestigate the same competition in the context of the FFTFIM. 
The main advantage is that, as we shall see, the FFTFIM has 
a plaquette phase in a simple limit, that of classical spins. A 
combination of large scale QMC simulations and semiclassi- 
cal arguments shows that the plaquette phase disappears when 
quantum fluctuations are introduced, and that columnar order 
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FIG. 1. (color online). Visualization of the square lattice structure 
(black solid lines for ferromagnetic couplings and dashed lines for 
AF couplings) and of particular classical spin configurations with 
corresponding dimer configurations on the dual lattice: (a) plaquette 
dimer structure and (b) a columnar dimer structure. Dimer densities 
(thick color lines) on the dual lattice correspond to bond energies on 
the original spin lattice as described in the text. 

is stabilized from the paramagnetic phase down to zero trans- 
verse field, hence also for the QDM for V/t = 0. 

A convenient gauge-invariant description of ordering in the 
FFTFIM can be introduced [16, 17] by using the energies of 
bonds {ii + e„) of the original lattice to define dimer densities 
di,a -5(1" J ii+e„ I J ) on the bonds of the dual lattice 

crossing them. This mapping gives rise to the visualization 
of two typical dimer-solid states shown in Fig. 1. It allows 
one to characterize all aforementioned dimer-solid phases by 
a single complex order parameter [11, 12, 15, 18, 19], which 
in terms of da{(\) - jj 2, the Fourier transforms of 

dimer density di^a, has a very simple form, 

^2[dA0,7r) + idy(7T,0)]. (3) 

Long-range ordering into a valence bond solid is indicated 
by (I*?!) > 0. However, the mere fact of a non-vanishing ex- 
pectation value for \^'\ tells us nothing about the precise struc- 
ture of the solid. To make such a distinction, one needs infor- 
mation about the phase of the order parameter = I'I'le'* 
[15]. In this language, the columnar phase is characterized 
by a purely real or imaginary 'I', i.e., (p - - nn/2 with 
n = 0, . . . , 3 and the plaquette state is described by a phase 
angle <p = (pp = jijA + nn/2 with n = 0, . . . , 3. A mixed-state, 
on the other hand, is an interpolation between these two sit- 
uations. It is a state that breaks translation and rotation sym- 
metry but where the phase (p is neither given by (p^ nor <pp. 
Instead of directly studying the phase (f>, we find it more con- 
venient to look at a phase sensitive order parameter [20, 2 1 ] 
defined, e.g., as 

^4 = <mcos(4.^)). (4) 

It is straightforward to see that in a columnar phase *I'4 \^'\ 
whereas *I'4 — > -\^'\ holds for plaquette ordering. In a true 
mixed phase, however, |*I'4| will not converge towards I*!*! in 
the thermodynamic limit. 

In addition, one may resort to special plaquette correlation 
functions P+ and P , as described in [14, 15], to distinguish 
the type of ordering. We will show however that in the present 
context those have to be analyzed with care. 



0.8 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
r 



FIG. 2. (color online). Plot of the order parameter -1*4 for the clas- 
sical limit of (1) at linear lattice size L = 16 and various low tem- 
peratures confirming the nature of the ground state to be a plaquette 
state. For comparison the norm yV] is shown at T = 0.005 as well 
as calculated from the exactly known ground state. The theoretical 
critical field Fc = 2 V2 is indicated by the vertical line. 

Classical limit. — The classical limit of Eq. ( 1 ) defined by 
interpreting cr~ and cr* as components of an ordinary three- 
dimensional vector of length S has previously been used to 
study the FFTFIM on the honeycomb lattice [17], which can 
be mapped onto the QDM on the triangular lattice in the limit 
of weak field, and its predictions, when combined with linear 
spin-wave theory (LSWT), agree very well with QMC results 
on the QDM. In particular, the phase adjacent to the paramag- 
netic phase is the 12-site plaquette phase revealed by QMC. 

For the FFTFIM on the square lattice, the paramagnetic 
phase with spins along the x direction can be shown to ex- 
tend down to Fc = 2 V2, and in the range < F < 2 V2, a 
single phase is stabilized with a four-site unit cell and one 
spin per unit cell fully polarized along the field. This has 
been checked numerically and can be proven analytically in 
the limits F — > and F — > F^ [22]. This phase coiTesponds 
to a plaquette state in the dimer language (see Fig. 1(a)). The 
columnar state would correspond to the dimer density pattern 
sketched in Fig. 1(b). Minimizing the energy for a unit cell 
which has the symmetry of the columnar state leads to a con- 
figuration which is a saddle point of the total energy, but not a 
local minimum, and a fortiori not a global one. 

This result is quite surprising in view of the most recent 
results obtained for the QDM at V/t = 0, which suggest a 
columnar or a mixed phase, but not a plaquette phase. In or- 
der to confirm it and to test the order parameter T introduced 
above, we have performed MC simulations (of the classical 
version of (1)) on systems of size L with periodic boundary 
conditions. They fully coiToborate our analysis, see Fig. 2. In 
particular, we find a convergence of the phase sensitive quan- 
tity -*I'4 — > l^l in the ordered region F < Fc, proving the 
presence of the plaquette state in the classical FFTFIM. For 
more details, see [22]. 

The stacked Ising model. — Motivated by this somehow 
unexpected result, and as a first step towards a full quantum 
treatment of the FFTFIM, we have decided to use the same 
type of MC simulations to revisit the stacked fully frustrated 
Ising model first studied by Jalabert and Sachdev [18]. This 
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FIG. 3. (color online), (a) Order parameter distribution function (in 
the complex T-plane) for the stacked magnet at T = 1.5 for vari- 
ous lattice sizes, (b) Finite-size behavior of the two order parame- 
ters and *I'4 versus 1/L. (c) Same as (b) for the structure factors 
P+ (plaquette correlations) and P_ (columnar correlations) defined 
in [14]. The color shaded regions indicate the approximate range 
where almost perfect l/(l)-symmetric distributions are found. All re- 
sults demonstrate the presence of columnar dimer solid in the stacked 
magnet. 
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FIG. 4. (color online), (a) Spin structure factor 5 (0, ;r) as a function 
of the transverse field F indicating an ordered phase below F < 1.6 
(Inset : Scaling of S(0,n) with the critical exponents from the XY 
universality class.) (b) Associated Binder parameter. It gives a pre- 
cise estimate of F^, indicated by the vertical line, (c) Scaling of the 
order parameter ratio ^4/^ for the quantum FFTFIM for various val- 
ues of the field F. It provides convincing evidence of columnar or- 
dering for the smallest fields, (d) Analysis of the structure factors P+ 
and P- used in Ref. [14] for the smallest field studied F = 0.01. 



model is closely related to the FFTFIM since they can be 
mapped onto each other using Trotter-Suzuki's decomposi- 
tion if appropriately scaled coupling constants and stack di- 
mension are used, and Jalabert and Sachdev have concluded 
that columnar order is stabilized. However, what they actually 
showed is that there is a phase transition at T ^ 2.85 below 
which 1^1 acquires a finite value without actually resolving 
the phase of the order parameter. We have performed high- 
statistics MC simulations inside the ordered phase at T = 1 .5, 
with a focus on obtaining such conclusive phase-sensitive in- 
formation. The resulting distribution functions of the complex 
order parameter for linear lattice sizes L - 16, 32, 64, 96 are 
shown in Fig. 3(a). Interestingly, a t/(l) symmetric distribu- 
tion is found for L < 32, but ultimately four peaks develop at 
phase angles = (pc and a clear columnar state emerges above 
the U(l) regime. These results lead to two important conclu- 
sions: (ii) Quantum fluctuations seem to change the nature of 
ordering from plaquette to columnar; (ii) A U(l) length scale 
(as in [21]) prohibits to draw any conclusion from data on 
small system sizes. 

To get a better understanding of the effect of this length 
scale, it is instructive to study its effect on the finite-size 
scaling of order parameters. Figure 3(b,c) compares the finite- 
size dependence of \^'\ and ^'4. It is only for L > 32 that a 
clear signal appears showing that will eventually converge 
towards \^'\. We have also studied the scaling of the two corre- 
lation functions P+ and P - introduced in Ref. [14] - meant 
to capture plaquette and columnar correlations separately (see 
Fig. 3(c)). Quite remarkably, the columnar correlations first 
go through a minimum upon increasing the size before ulti- 
mately reaching a finite-expectation value for L oa. By 



contrast, the plaquette correlations seem to tend to a finite ex- 
pectation value before ultimately dropping abruptly to zero. 
These results show that the U{1) regime has a strong influ- 
ence on finite-size scaling, and that *I'4 is probably more reli- 
able than P+ and P : For small sizes, the sign of *i'4 is already 
correct, while the tendencies for P+ and P- are completely 
misleading. 

The quantum FFTFIM. — Next, we turn to an exact treat- 
ment and discuss the ground-state phase diagram of Eq. ( 1 ) 
for arbitrary field F. Using truly unbiased quantum Monte 
Carlo methods [23, 24] on periodic lattices with linear sys- 
tem sizes up to L = 256, we investigate the order parameter 

and the basic ordering modes in terms of the (sublattice) 
structure factors 5(0,0) and S{Q,n) (see. Fig. 1(b)). In each 
case, the inverse temperature fi was carefully chosen to be 
large enough to sample ground-state properties. Results for 
the spin structure factors are reproduced in Fig. 4(a), which 
show that an ordered phase is stabilized by quantum fluctu- 
ations below F < 1.6. A more precise estimate of the crit- 
ical transverse field is obtained from the Binder parameter 
g2 = 1 - <5(0,7r)2)/(3(5(0,7r))^) which is scale invariant at 
the quantum critical point. The crossings of Q2 lead to the 
estimate F^, = 1.578(2), see Fig. 4(b). This critical field is 
considerably reduced as compared to the classical value 2 V2. 
Using the value obtained for and 3D XY critical exponents 
(j6 and v) [25], a good data collapse of the structure factors is 
obtained, directly demonstrating the 3D XY universality class 
of the quantum critical point. This result coherently adds to 
the indirect analysis performed in Ref. [20] and supports the 
field-theoretic picture [26, 27] beyond reasonable doubt. It 
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also confirms the similarity between the triangular and square 
lattice FFTFIM as far as the critical behavior is concerned. 
An analysis of the FFTFIM at finite temperature that reveals 
the presence of an extended critical phase will be presented 
elsewhere [22]. To determine the nature of the phase(s) for 
r < Fc, we calculate the order parameter *i'4 and analyze its 
finite-size scaling with respect to Figure 4(c) shows the 
essential results of this analysis for a couple of selected field 
values inside the ordered region. While it is clearly difficult to 
get phase sensitive information close to the quantum critical 
point (F = 1 .2) - continuous distributions as in Fig. 3(top) are 
seen up to the largest system size - we note that it actually 
becomes easier to obtain structural information in the limit 
F/y — > 0. Indeed, for F - 0.01 the scaling of *I'4 is sufficient 
to conclude on the presence of a columnar state, a result fully 
supported by the behavior of P+ and P , see Fig. 4(d). Based 
on this analysis, we can safely conclude that the columnar 
state is stabilized at least up to F a; 1 .0. Above this value, it is 
currently difficult to analyse the nature of the phase by means 
of finite-lattice simulations. This particular difficulty close to 
the quantum critical point is actually not surprising since the 
3D XY quantum critical point has a C/(l) symmetry which is 
at the origin of the disturbing crossover behavior found here 
[21]. A similar observation has already been made in the case 
of the triangular lattice [20]. 

Semiclassical approach. — All these results indicate that 
the order in the low field phase changes from plaquette to 
columnar upon going from the classical to the ultra-quantum 
spin-1/2 case. In order to investigate this issue from a com- 
plementary angle, we have studied how quantum fluctuations 
may favor columnar order in the context of a semiclassical 
approach. As stated earlier, the classical solution which cor- 
responds to columnar order is not a local minimum of the en- 
ergy. So the spectrum of the harmonic fluctuations cannot 
be diagonalized by a Bogoliubov transformation. An upper 
bound of the zero point energy can nevertheless be obtained 
by applying local fields in the direction of the spins to achieve 
the positiveness of the quadratic form that can be diagonalized 
by a Bogoliubov transformation (see Supplemental Material). 
For the plaquette state, the Hamiltonian of the quadratic fluc- 
tuations can be diagonalized by a Bogoliubov transformation. 
The phase diagram obtained by comparing the energy of the 
columnar and plaquette states corrected in this way by quan- 
tum fluctuations is shown in Fig. 5. The results nicely confirm 
the anticipated transition: upon reducing the spin, quantum 
fluctuations reduce the energy of the columnar much faster 
than that of the plaquette state and stabilize it for all fields for 
spin-1/2. 

Conclusions. — From the various results reported in this 
paper, a number of firmly established conclusions can be 
drawn. First of all, at low field, the ground state of the quan- 
tum FFTFIM on the square lattice is an ordered phase which, 
in the gauge invariant dimer language, corresponds to colum- 
nar order This order is stable at least up to F/y - 1, and 
possibly up to the transition to the paramagnetic phase, which 
has been shown to take place at TjJ ^ 1.58. By contrast. 
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FIG. 5. (color online). Phase diagram in the F - 1/5 plane, showing 
the phase boundary between columnar and plaquette phase where 
the transition line shown is an upper bound for the stability of the 
plaquette phase. The columnar phase is clearly stabilized by quantum 
fluctuations. 



we have shown that the classical version of the model de- 
velops an order with a four-spin unit cell at low temperature 
which, in the dimer language, corresponds to plaquette order. 
Furthermore, the robustness of plaquette order with respect 
to two types of perturbations, stacking and quantum fluctua- 
tions, has been tested, with the conclusion that it is relatively 
fragile: 3D stacking clearly stabilizes columnar order, while 
quantum fluctuations treated within linear spin-wave theory 
stabilize columnar order in a large portion of the (1/5, F/y) 
phase diagram, and in particular for all values of F/7 outside 
the paramagnetic phase for the quantum model with S = 1/2. 

Regarding the closely related QDM on the square lattice, 
the present results clearly demonstrate that, in the purely ki- 
netic regime V/f = 0, which is equivalent to the F/7 ^ limit 
of the FFTFIM, columnar order is stabilized. In addition, this 
is likely to be the case in a finite parameter range around this 
point since we have not seen any particular indication of a 
proximity to a critical point. Complementary work carried 
out in parallel to ours comes to similar conclusions [28]. 

Finally, the other important message of this paper is the 
presence of a surprisingly large length scale below which the 
model has an effective U(l) symmetry. This makes the iden- 
tification of the actual order stabilized in the thermodynamic 
limit particularly difficult and, to a large extent, explains the 
conflicting results obtained so far on the type of order stabi- 
lized in the QDM between V/f - and Vjt - 1 . By study- 
ing very large systems, we have been able to beat this length 
scale up to F/y ^ 1. Beyond that value, the proximity of the 
transition to the paramagnetic phase makes this length scale 
even larger, and whether columnar order remains stable up to 
the paramagnetic phase, as suggested by the semiclassical ap- 
proximation, could not be numerically decided. 
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Supplemental Material 

Here we present the details of the analysis of the effect of 
quantum fluctuations for the plaquette and columnar states. 
Quantum fluctuations are investigated in the context of a large 
S expansion (spin wave expansion). In this framework, the 
fluctuations around a classical solution are described in terms 



of Holstein-PrimakofF bosons and the fluctuation Hamiltonian 
is truncated at the harmonic level (linear spin wave expan- 
sion). 

In the two classical structures that are compared, the fluc- 
tuation Hamiltonians do not contain terms which are linear 
in the Holstein-Primakoff bosons and therefore are purely 
quadratic. In the case of the plaquette state, the correction to 
the classical energy is obtained by diagonalizing the fluctua- 
tion Hamiltonian. This is done in Fourier space via a standard 
Bogoliubov transformation which yields the dispersion rela- 
tions and the correction to the classical energy. 

The computation in the case of the columnar state is more 
involved. This time the starting classical configuration is not 
a minimum but a saddle point of the classical energy. This re- 
sults in a quadratic fluctuation Hamiltonian which is not posi- 
tive definite. Trying to diagonalize the quadratic Hamiltonian 
via a Bogoliubov transformation in this case yields a spectrum 
which is not well defined for all values of momenta. For the 
regions in the Brillouin zone near k = (0, ±7r) the spectrum 
is not real but takes complex nonphysical values. In this case, 
the correction to the classical energy cannot be computed. The 
fact that the spectrum is not well defined results from the trun- 
cation of the spin wave approximation to harmonic level. If 
this state is to become the ground state when quantum fluc- 
tuations are included, higher order terms in the spin wave ex- 
pansion must yield a spectrum with real frequencies. However 
this implies treating higher order terms self consistently. 

In order to have a well-defined spectrum we proceed in a 
different way. We add to the Hamiltonian a local field term on 
each site of the form 

j 

where 6 > parametrizes the local field strength and where 
the summation is taken over all sites and at each site the axis 
z'ij) is directed along the direction of the classical spin. The 
field 6 is adjusted to have a fluctuation Hamiltonian which 
is positive definite allowing it to be diagonalized by a Bo- 
goliubov transformation. The resulting spectrum has real and 
positive frequencies with soft modes only at the wave vectors 
k - (0, ±n). The advantages of this approach are the follow- 
ing: the addition of V to the Hamiltonian does not change the 
classical energy of the state considered and allows to obtain 
dispersion relations which are physically meaningful. Fur- 
thermore y is a strictly positive contribution to the Hamil- 
tonian. Hence the corrections to the energy of the columnar 
state computed with this approach provide an upper bound of 
the energy of this state at order 1/5 . 

We find that as a function of S and TjJ, the upper bound 
of the energy of the columnar state is lower than the energy of 
the plaquette state in a significant parameter range (see Fig. 5 
of main text). In particular, the columnar state wins for all 
fields for 5 = 1/2, the value of the spin for which the model 
can be exactly mapped onto a QDM with V/f = in the limit 

r/y ^ 0. 



